Compositional (version 4.3)

Contour plot of a Flexible Dirichlet distribution in S^2: Contour plot of a Flexible Dirichlet distribution in \(S^2\)

Description

Contour plot of a Flexible Dirichlet distribution in \(S^2\).

Usage

fd.contour(alpha, prob, tau, n = 100, x = NULL)

Arguments

alpha

A vector of the non-negative \(alpha\) parameters.

prob

A vector of the clusters' probabilities. It must sum to one.

tau

The non-negative scalar \(tau\) parameter.

n

The number of grid points to consider over which the density is calculated.

x

This is either NULL (no data) or contains a 3 column matrix with compositional data.

Value

A ternary diagram with the points and the Flexible Dirichlet contour lines.

Details

The user can plot only the contour lines of a Dirichlet with a given vector of parameters, or can also add the relevant data should they wish to.

References

Ongaro A. and Migliorati S. (2013). A generalization of the Dirichlet distribution. Journal of Multivariate Analysis, 114, 412--426.

Migliorati S., Ongaro A. and Monti G. S. (2017). A structured Dirichlet mixture model for compositional data: inferential and applicative issues. Statistics and Computing, 27, 963--983.

See Also

norm.contour, bivt.contour, comp.kerncontour, mixnorm.contour

Examples

# NOT RUN {
x <- rfd(n = 50, a = c(12, 11, 10), p = c(0.25, 0.25, 0.5), tau = 4 )
mod <- fd.est(x, ini.iter = 10, final.iter = 20)
fd.contour(mod$alpha, mod$prob, mod$tau, x = x)
# }